Heavy-quark correlations in deep inelastic scattering 
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We discuss some results for heavy quark correlations in next-to-leading order in deep inelastic electroproduction. 



1. INTRODUCTION 

Order a s QCD corrections to structure func- 
tions containing heavy quarks and to single heavy 
quark inclusive distributions in deep-inelastic 
electroproduction {i.e., r J*{q) + P{p) — > Q(pi) + 
X where X stands for any final hadronic state al- 
lowed by quantum- number conservation and P(p) 
is a proton of momentum p ) were recently pub- 
lished in jl) and respectively. By combin- 
ing the next-to-leading order (NLO) heavy quark 
structure functions with the corresponding light- 
quark structure functions j|, the heavy quark 
content of the nucleon has been studied as a func- 
tion of Q 2 = -q 2 and x = Q 2 /2p ■ q j§. Event 
rates for charm production integrated over bins 
in x and Q 2 relevant to HERA data have been 
calculated in ||. 

To further the study of deep-inelastic electro- 
production of heavy quarks we have recalculated 
the virtual-photon-parton cross sections of in 
an exclusive fashion || . This enables us to study 
the single and double differential distributions 
and correlations among all outgoing particles in 
the virtual-photon induced reaction 7* + P — > 
Q + Q + X with X = or 1 jet and to easily incor- 
porate experimental cuts. By examining distri- 
butions and correlations that are trivial at lowest 
order (for example, the azimuthal angle between 
the heavy quark and the heavy antiquark) one di- 
rectly tests the hard scattering cross section that 
is predicted by perturbative QCD and factoriza- 
tion. Therefore our results should give a clean 
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test of NLO perturbative QCD. Nevertheless we 
remind the reader that this is a fixed order per- 
turbative calculation and suffers from the same 
problems as all NLO calculations. There are re- 
gions in phase space where it will break down. 
For example, in the above mentioned azimuthal 
angle distribution, if one looks at the prediction 
too near the back-to-back configuration there will 
be an extra enhancement of the cross section due 
to multiple soft gluon emission which is not in- 
cluded in our NLO result. 

In this short report we discuss some details of 
the calculation and present interesting distribu- 
tions. The transverse and longitudinal photon 
components are treated separately and the lat- 
est CTEQ3 parton densities 0, consistent with 
the newly released HERA data ||, are used in 
the kinematic regime appropriate for production 
of charm quarks at HERA. We make our predic- 
tions at fixed values of Q 2 (> 8.5(GeV/c) 2 ) and 
x (> 4.2 x 10~ 4 ). 

We stress that here we only consider extrinsic 
heavy quark production, involving Bethe-Heitler 
and Compton production from massless partons. 
For a discussion of intrinsic production, where the 
heavy quark is considered to be part of the pro- 
ton's wavefunction, see Brodsky et al. || . A vari- 
able flavor scheme which joins the extrinsic heavy 
flavor production picture at /iphys *C m with a 
light mass parton density picture at /i P hvs ^ m 
has been discussed by Aivazis et al. (n]] . By com- 
paring the fixed flavor scheme calculation of W 
with the variable flavor number scheme of |l0[ 
it is concluded in that the former yields the 
most stable and reliable results for £2(2;, Q 2 , m 2 ) 



in the threshold region ( i.e., Q 2 < 10m 2 where 
m is the mass of the heavy quark) . The contribu- 
tion by W.L. van Neerven to these Proceedings 
reviews the progress towards implementing the 
variable flavour number scheme in NLO 

In our exclusive computation we use the sub- 
traction method which is based on the replace- 
ment of divergent (collinear or soft) terms in the 
squared matrix elements by generalized plus dis- 
tributions. This allows us to isolate the soft and 
collinear poles within the framework of dimen- 
sional regularization without calculating all the 
phase space integrals in a spacetime dimension 
n ^ 4 as usually required in a traditional single 
particle inclusive computation. In this method 
the expressions for the squared matrix elements 
in the collinear limit appear in a factorized form, 
where poles in n — 4 multiply splitting functions 
and lower order squared matrix elements. The 
cancellation of collinear singularities is then per- 
formed using the factorization theorem jl3| . The 
expressions for the squared matrix elements in 
the soft limit appear in a factorized form where 
poles in n~ 4 multiply lower order squared matrix 
elements. The cancellation of soft singularities 
takes place upon adding the contributions from 
the renormalized virtual diagrams. Since the final 
result is in four-dimensional space time, we can 
compute all relevant phase space integrations us- 
ing standard Monte Carlo integration techniques 
|l4| ] and produce histograms for exclusive, semi- 
inclusive, or inclusive quantities related to any 
of the outgoing particles. We can also incorpo- 
rate any reasonable set of experimental cuts. Our 
computer code has no small phase space slice pa- 
rameters and the parameters defining the gener- 
alized plus distributions may be tuned to give fast 
numerical convergence, which is an advantage of 
using this subtraction method. 

2. The 7*g channel 

Analysis of the partonic reaction 

7*(<z) + q(h) -> q(k 2 ) + Q(pi) + Q(P2) , (1) 

does not involve soft or virtual contributions so 
we can use it to most simply explain the method. 
One particular graph out of the order eg 2 set is 



Figure 1. One order eg 2 diagram contributing to 
the amplitude for the reaction (1). Additional 
graphs are obtained by reversing the arrows on 
the light-quark lines (dashed). 



shown in fig. 1 . We project this set of graphs 
on a particular photon polarization state i and 
square the amplitude. The answer is written as 
Cj ig M?(3) where the coefficient absorbs overall 
factors. Now we write the photon-parton cross 
section as 



M% = C i<q M?(3)dT 3 



(2) 



where the 3 indicates a three-body final state. To 
proceed we need some notation for the variables 
s = 2ki -q,t\ = t — m 2 = —2k\ -p 2 , u\ — u\~q 2 = 
2q ■ p2, t' = —2ki ■ k 2 , U 1 — q 2 — 2q ■ k 2 and Sg = 
(pi + P2) 2 — q 2 ■ Then if we introduce x = s^/s 
and cos -1 y — q ■ k 2 j\q ■ k%\ we can write t' = 
-(l/2)s'(s'/s)(l -x)(l + y). This representation 
shows that t' — > as y — > — 1. This t' divergence 
yields the collinear divergence in M?(3), which 
is due to the presence of the Feynman graph in 
fig.l. To handle this divergence we multiply and 
divide by t' so that 



da® = C i<g f?(x,y,6 u 6 2 )dr 3 , 
where 



(3) 



f?(x,y,6 l: e 2 )=t'M?{Z), (4) 
is now finite as y — > — 1 and 

dr; = dT 3 /t'. (5) 

is divergent. We next replace the factor 
(l+J/r 1+e/2 = (l/(l+y)) w +*(l+tf)(l/e+lnw),(6) 



in the (n-dimensional) three-body phase space, 
where < u> < 2. This has the effect of splitting 
the n-dimensional cross section into two divergent 
parts. However it lets us take the limit e — > and 
get back to 4-dimensions, with u> regularizing the 
collinear singularity. Replacing the divergent fac- 
tor (\ + y)- 1+t ' 2 one obtains the following 
decomposition: (for more details see M) 



with 
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and 

daft = --C hq HNdT 2 ( S r /2 (-Y /2 

7T \ S / 

V 1 + y/ui 

x dy sin £ 9 2 d9 2 f?(x,y, 9 U 2 ). (9) 
The function / satisfies 

f?( x , -l, Oi,e a ) = f!(x, Ox) + f?(x, di,9 2 ) , (io) 

with 

(11) 



f?(x,9 u 9 2 ) sin 6 9 2 dB 2 =0. 
For reaction (1) we find that 
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where Bi^qed it the Born amplitude with k\ re- 
placed by xfci. The factor in the square brack- 
ets is the Altarelli-Parisi splitting function in n 
dimensions so we can define new (collinear free) 
cross sections 
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l, 02 ) are available 



As the quark channel only contains collinear 
singularities, we also use it to illustrate how the 
generalized plus distributions are implemented 
numerically and how the u) dependence disap- 
pears in the sum (pit). To this end consider 
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where we have suppressed all indices and argu- 
ments of ff(x, y, 9i,9 2 ) other than y. In the bot- 
tom line we see that the infinity encountered at 
the lower integration limit y = — 1 is cancelled in 
the sum of the two integrals. In practice one can 
only reasonably take the lower limit to be —1 + S 
where S ~ 10~ 7 in double precision Fortran be- 
fore round off errors enter. None the less, the 
final results are stable with respect to the varia- 
tion of 8 in the range 10~ 5 to 10~ 7 . The upper 
limit of the second integral gives a contribution 
/(— 1) Inw which cancels against the In a; term in 
dof q . The first integral in the bottom line is 
commonly called the "event" and has a positive 
definite weight. The second integral plus the fac- 
torized collinear contribution is commonly called 



the "counter-event" and may have either positive 
or negative weight. 

The implementation of this procedure for the 
gluon channel is similar but more complicated due 
to the presence of both soft and collinear diver- 
gences. 



3. Results 

Recalling that the probability density is re- 
lated to the momentum density via fi/ni^i M/) = 
M/) we write the hadronic cross section 

as 



(17) 



We now specialize to the case where H is a proton, 
as in the case of HERA. Using the relations 



F, 



Q 2 

An 2 a 



(18) 



where k = 2, L with er 2 = ctg + 3ctl/2, and rela- 
tions for the scaling functions, we find 



F k (x,Q 2 ,m 2 



Q 2 a s (p 2 
47r 2 m 2 
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c (1) ln-^ 
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where k = 2, L. We have set /j,/ = /i and shown 
the /i 2 dependence of the running coupling a s ex- 
plicitly. The lower boundary on the integration is 
given by £ m ; n = ie(4to 2 + Q 2 )/Q 2 . This formula 
yields the standard heavy quark structure func- 
tions F 2 (x,Q 2 ,m 2 ) and Fl(x, Q 2 : m 2 ) for elec- 
tron proton scattering, and we will present results 
as differentials of these functions. From the for- 
malism described in the previous section we are 
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Figure 2. The distributions dF 2 (P t )/dP t for 
charm- ant icharm pair production at fixed Q 2 = 
12(GeV/c) 2 with x = 4.2 x 1CT 4 (solid line), 
8.5 x 1CT 4 (dotted line), 1.6 x 1(T 3 (short dashed 
line) and 2.7 x 10 -3 (long dashed line). 



left with events of positive weight and counter- 
events of either positive or negative weight. Our 
program outputs the final state four vectors of the 
event and counter-event together with the corre- 
sponding weight. We histogram these into bins 
to produce differential distributions. 

We start by considering the production of 
charm quarks. We use m = m c = 1.5GeV/c 2 
and simply choose the factorization (renormaliza- 
tion) scale as fi 2 = Q 2 + 4(m 2 + (P t c + P t E ) 2 /4). 
Note that there are many possible choices of scale 
as we have all components of the final four vec- 
tors. Aside from the Pt dependence, this choice 
reduces to /i 2 = Q 2 for electroproduction of mass- 
less quarks and /i 2 = 4to 2 for the photoproduc- 
tion of charm quarks. We introduce a Pt depen- 
dence by adding in the average of the magnitude 
of the transverse momenta of the heavy quark and 
heavy antiquark. As mentioned earlier we use the 
CTEQ3M parton densities @ in the MS scheme 
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Figure 3. The distributions dF 2 (P t )/dP t for 
charm-anticharm pair production at fixed x = 
8.5 x 10" 4 and Q 2 = 8.5 (solid line), 12 (dot- 
ted line), 25 (short dashed line), 50 (long dashed 
line) all in units of (GeV/c) 2 . 



and the two loop a s with A4 = 0.239 GcV. 

The first distribution we present basically mea- 
sures the transverse momentum of the addi- 
tional jet which recoils against the heavy quark 
pair. The Pt distribution of the charm-anti- 
charm pair is shown in fig. 2 where we 
plot dF 2 (x,Q 2 ,m 2 ,Pt)/dP t as a function of P t . 
The histograms are presented at fixed Q 2 — 12 
(GeV/c) 2 for x values of 4.2 x 10~ 4 , 8.5 x 10" 4 , 
1.6 x 10~ 3 and 2.7 x 10~ 3 respectively. One sees 
that the Pt distribution peaks at small Pt and has 
a negative contribution in the lowest bin. This 
is a region where the dominant contribution is 
from counter-events so the weights can be nega- 
tive. The results of this calculation require miss- 
ing contributions from even higher order pertur- 
bation theory before this bin will have a posi- 
tive weight. The depth of the negative bins is a 
function of x, Q 2 , and the choice of scale. Note 
that at larger P t the structure function is domi- 



nated by the contribution from the square of the 
bremsstrahlung graphs so the weights are posi- 
tive. Figure 3 shows the corresponding results 
for fixed x = 8.5 x 10~ 4 plotted for the Q 2 val- 
ues of 8.5 (GeV/c) 2 , 12 (GeV/c) 2 , 25 (GeV/c) 2 
and 50 (GeV/c) 2 respectively. The distributions 
peak near small Pt and are either small or nega- 
tive in the first bin. The histograms with Q 2 = 
12(GeV/c) 2 and i = 8.5x 10" 4 (the dotted line) 
are the same in figs. 2 and 3. We have also used 
the same scales on the axes so one can easily see 
that there is a greater change if we fix x and vary 
Q 2 than if we fix Q 2 and vary x. We will continue 
to use the same scale for all the pairs of later plots 
to simplify the comparison between them. 

We now turn to the distributions in the az- 
imuthal angle between the outgoing charm quark 
and charm antiquark which we denote as Aip. 
This is the angle between the P< vectors of the 
heavy quark-antiquark in the c. m. frame of the 
virtual-photon-hadron system. Since we integrate 
over the azimuthal angle between the plane con- 
taining the incoming and outgoing leptons and 
the plane defined by the incoming parton and 
outgoing heavy quark (to limit our discussion to 
F2 and Fl) we can only plot relative azimuthal 
correlations. In the Born approximation this dis- 
tribution is a delta function at ir, as their four 
momenta must balance. Due to the radiation of 
the additional light mass parton, the distribution 
has a tail extending below it and has a valley at n. 
The distributions become negative in the highest 
bins. This negative region is a general feature of 
all exclusive calculations. Figure 4 shows results 
for dF 2 (x, Q 2 , m 2 c ,A(j))/d(A(j>) at the same values 
of fixed Q 2 and variable x as chosen previously 
in fig. 2, while fig. 5 shows the results for fixed 
x and variable Q 2 as chosen previously in fig. 3. 
Note again that the dotted histograms are the 
same in figs. 4 and 5, and there is more variation 
for fixed x and changing Q 2 than for fixed Q 2 and 
changing x. 

In this report we have outlined the NLO calcu- 
lation of the virtual-photon-parton (Wilson) co- 
efficient functions in the fully differential produc- 
tion of heavy quarks plus one jet. More correla- 
tions and checks are given in |@|. The computer 
program we have written for the exclusive cal- 
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Figure 4. The distributions dF 2 {A(j)) / d(A</>) for 
charm-anticharm pair production at the x and Q 2 
values given in Figure 2. 



culation has the advantage that the four vectors 
of the heavy quark, heavy- antiquark and/or one 
additional light parton jet are produced for each 
event and can be subjected to experimental cuts. 
We have not done this in any of the plots shown 
here but the computer program is available and 
can be easily modified to incorporate acceptances 
of the detectors at HERA. The program has been 
used to reanalyse the EMC data on charm pro- 
duction in |L7j . 

4. Conclusion 

In this report we have presented some plots 
which depend on information from the four vector 
of the additional jet. We showed the distributions 
in the transverse- momentum (P t ) of the heavy 
quark antiquark pair and in the azimuthal angle 
(A0) between the Pt vectors of the heavy quark 
and heavy antiquark. All quantities were pre- 
dicted in the c. m. frame of the photon- proton sys- 
tem after integration over the azimuthal angle be- 
tween the plane containing the incoming and out- 



12 3 

A* 

Figure 5. The same distributions at the x and Q 2 
values given in Figure 3. 



going lepton and the plane containing the incom- 
ing proton and outgoing heavy quark. The results 
were presented as distributions in F2(x,Q 2 ,m 2 ) 
at specific points in x, Q 2 and m 2 = m 2 . None of 
these distributions can be reproduced by any K- 
factor multiplication as the corresponding Born 
distributions are proportional to delta-functions. 
In all cases the histograms of these distributions 
have negative bins. These are regions where the 
NLO calculation is not sufficient and a NNLO or- 
der calculation (or some form of resummation) 
should be made. A general statement about the 
magnitude of the NLO contribution compared to 
the LO one is difficult to make as the size and sign 
of the corrections may vary strongly between dif- 
ferent regions of phase space. However, we see 
that all plots have larger Q 2 variation at fixed 
x as comparied to varying x at fixed Q 2 . By 
varying the renormalization (factorization) scale 
we observed that the distributions presented here 
changed in normalization but not in shape. In our 
study we also calculated averages of various quan- 
tities and found the typical variation between cen- 



tral and extreme scale choices of 5 percent for 
charm production. 
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